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ABSTRACT 


The  performance  of  routine  FNWC  sea  level  pressure  anal- 
yses and  forecasts  in  the  Beaufort  Sea  is  evaluated  for  driv- 
ing an  ice  model.  The  geostrophic  winds  determined  from  the 
FNWC  pressure  fields  are  compared  to  geostrophic  winds  de- 
termined from  the  AIDJEX  experiment  observations.  A simu- 
lation with  the  AIDJEX  model  has  been  run  using  FNWC  sea  level 
pressure  analysis  fields  for  the  period  15  to  25  May  1975 
and  compared  to  the  AIDJEX  results. 

RMS  sea  level  pressure  errors  of  about  3 mb  were  found 
with  analysis  errors  reduced  by  about  1 mb  when  one  AIDJEX 
observation  was  included  in  the  analysis.  Geostrophic  winds 
had  speed  differences  of  about  40%  and  direction  differences 
of  about  50°.  The  FNWC  simulation  with  the  AIDJEX  model 
produced  a position  error  of  about  7 km  for  one  station 
over  the  ten  day  run. 

Recommendations  are  made  to  improve  the  FNWC  hemispheric 
fields  in  the  Arctic  for  application  to  ice  forecasting. 
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I.  INTRODUCTION 

Recent  developments  have  demonstrated  the  strategic 
importance  of  the  Arctic  to  the  security  of  America.  Soviet 
Delta-class  submarines  are  on  patrol  out  of  their  homeport 
of  Murmansk  armed  with  nuclear  tipped  missiles  able  to  strike 
the  heartland  of  America  from  anywhere  within  the  Arctic  Ocean. 

Soviet  strategy  is  to  seek  open  spots  in  the  ice  pack  as 
launch  points  for  their  attacks  [Anonymous,  1977].  The 
Northern  Sea  Route  represents  a major  Soviet  strategic  inter- 
est [Synhorst,  1973]  . As  reported  in  the  Monterey  Peninsula 
Herald,  (1977) , the  Soviets  have  demonstrated  their  ability 
to  penetrate  to  the  North  Pole  with  an  atomic  ice  breaker  and 
announced  their  intention  to  open  a great  circle  shipping 
route  through  the  Arctic  Ocean,  events  which  make  their  inter- 
est all  the  more  ominous.  German  submarines  in  World  War  II 
demonstrated  a threat  to  the  Northern  Sea  Route  by  their 
operations  in  the  Kara  Sea  [Synhorst,  1973].  In  the  event 
of  war  any  attempts  by  U.  S.  submarines  to  counter  the  Soviet 
SSBN  threat  in  the  Arctic  or  sever  the  Soviet  lines  of  com- 
munications via  the  Northern  Sea  Route  must  be  aided  by  re- 
liable knowledge  of  ice  conditions  in  order  to  make  acoustic 
performance  predictions  and  optimum  under-ice  track  routing. 

The  North  Slope  of  Alaska  holds  a critical  supply  of 
America's  lifeblood,  oil.  The  Navy  provides  forecasts  of  ice 
conditions  vital  to  the  continued  free  passage  of  U.  S. 
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maritime  resupply  vessels  to  this  valuable  area.  Again,  re- 
liable ice  forecasts  are  crucial  to  America's  interests. 

A.  NATURE  OF  THE  PROBLEM 

A knowledge  of  ridge  and  lead  formation  and  distribution 
are  of  importance  in  meeting  several  Navy  requirements.  For 
example,  the  long,  smooth  re frozen  leads  are  the  basic  landing 
strips  on  the  pack.  Knowing  the  location  and  ice  thickness 
of  these  refrozen  leads  would  be  of  great  help  in  carrying 
out  scientific  research  and  logistic  support  on  the  pack.  In 
shallow  areas  such  as  the  Chukchi  Sea  the  ice  keels  beneath 
ridges  can  form  a considerable  hazard  to  under-ice  navigation. 
Ice  free  leads  or  thin  ice  can  act  as  places  for  emergency 
surfacing  of  a submarine.  Penetration  of  the  ice  is  impeded 
when  the  pack  is  under  high  internal  stress.  Combined  with 
a knowledge  of  ice  thickness,  a knowledge  of  the  stress  state 
of  the  ice  would  aid  in  these  decisions,  and  weapons  and 
sensor  deployment  through  the  pack.  A submariner  might  find 
accurate  forecasts  of  these  conditions  very  reassuring  on  a 
submerged  transit  of  the  Arctic. 

The  two  features  in  the  polar  environment  that  most 
strongly  influence  underwater  sound  are  the  permanent  ice 
cover  and  the  sound  speed  structure  of  the  water.  The  essen- 
tially positive  sound  speed  gradient  causes  sound  rays  to  be 
generally  refracted  upward  and  then  reflected  and/or  scattered 
from  the  rough  ice-water  interface  [Diachok,  1976].  Both 
high  and  low  frequencies  are  rapidly  attenuated,  the  former 
by  scattering  losses  from  the  ice  cover,  the  latter  by  the 
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fact  that  very  low  frequencies  are  not  effectively  trapped 
in  the  channel  [Urick,  1975].  Recently  it  has  been  demon- 
strated that  back-scattering  strength  is  controlled  virtually 
exclusively  by  the  size  and  number  of  sea-ice  ridges  [Diachok 
1976]  . 

The  principal  source  of  ambient  noise  in  the  Arctic  and 
sub-Arctic  water  is  large  scale  deformation  of  ice  floes 
[Diachok,  1976].  When  the  ice  is  not  continuous,  as  near  the 
ice  edge,  noise  levels  5 to  10  dB  higher  than  those  measured 
at  the  same  sea  state  in  ice-free  water  have  been  observed 
[Urick,  1975].  The  mechanisms  responsible  are  related  to 
wind-generated  ocean  wave  interactions  with  sea  ice  at  the 
boundary  [Diachok,  1976].  Under  a continuous  ice  sheet  with 
calm  winds  and  rising  temperatures  very  low  levels  of  ambient 
noise,  well  below  the  lowest  Knudsen  curve  for  sea-state  zero 
have  been  found.  On  the  other  hand,  levels  some  40  dB  higher 
have  been  observed  when  the  ice  had  been  cracking  under  fall- 
ing air  temperature  [Urick,  1975] . The  acoustic  propagation 
and  ambient  noise  factors  are  both  important  in  under-ice 
navigation,  submarine  tracking  and  localization,  and  weapons 
performance. 

In  addition  to  its  effect  on  the  under-ice  environment 
the  ice  cover  has  a significant  effect  on  the  climate  of  the 
Arctic  by  acting  as  a lid  on  the  ocean.  A typical  mesoscale 
region  in  the  ice  pack  can  be  expected  to  contain  open  water, 
young  ice  tens  of  centimeters  in  thickness,  perennial  and 
deformed  ice  a few  meters  thick.  Of  concern  in  atmospheric 
modeling  is  the  fact  that  open  water  and  thin  ice  can  exert 


example,  results  by  Maykut  [1976]  indicate  ice  0-40  cm  thic 
allows  a net  heat  exchange  to  the  atmosphere  one  to  two  ord 
of  magnitude  greater  than  perennial  ice  a few  meters  thick. 
Ice  thickness  greater  than  100  cm  exhibit  little  change  in 
heat  flux  as  the  ice  thickens.  In  the  Central  Arctic  durin 
winter,  estimates  of  open  water  and  young  thin  ice  are  less 
than  1%  and  8 to  12%,  respectively.  Thus,  even  at  this  sma 
percentage  the  effects  of  these  two  ice  conditions  can  begi 
to  dominate  the  large  scale  heat  exchange  and  alter  boundar 
conditions.  Atmospheric  models  which  assume  the  ice  cover 
is  a continuous  three  meter  layer  may  seriously  underestimd 
the  heat  input  to  the  boundary  layer. 

At  the  present  time  the  Navy  ice  prediction  models  arai 
only  intended  to  serve  as  an  aid  to  hand  analysis.  The  wii 
drift  model  used  by  the  Navy  ice  forecasting  group  at  the 
Fleet  Weather  Facility,  Suitland,  Maryland,  gives  ice  spee< 
as  a linear  function  of  surface  wind  speed  [Skiles,  1968]. 
Moving  to  the  right  of  the  wind,  the  drift  angle  is  given  l 
an  inverse  exponential  function  of  wind  speed.  Under  de- 
velopment is  a model  to  provide  daily  wind  drift  after  .j 
Shuleikin  [1953]  and  growth  after  Zubov  [1945]  at  62  coast 
weather  reporting  stations  [Gerson  and  Simpson,  1976,  and 
Gerson,  1975],  Barnett  [1976]  has  presented  a scheme  for 
estimating  the  general  ice  conditions  during  the  shipping 
season  of  the  North  Slope  sea  route  using  the  height  of  t] 
1000  mb  surface  at  two  geographic  locations  (52N,  100E,  a; 
70N,  140E) . 
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a large  influence  on  regional  heat  exchange  rates.  For 
example,  results  by  Maykut  [1976]  indicate  ice  0-40  cm  thick 
allows  a net  heat  exchange  to  the  atmosphere  one  to  two  orders 
of  magnitude  greater  than  perennial  ice  a few  meters  thick. 

Ice  thickness  greater  than  100  cm  exhibit  little  change  in 
heat  flux  as  the  ice  thickens.  In  the  Central  Arctic  during 
winter,  estimates  of  open  water  and  young  thin  ice  are  less 
than  1%  and  8 to  12%,  respectively.  Thus,  even  at  this  small 
percentage  the  effects  of  these  two  ice  conditions  can  begin 
to  dominate  the  large  scale  heat  exchange  and  alter  boundary 
conditions.  Atmospheric  models  which  assume  the  ice  cover 
is  a continuous  three  meter  layer  may  seriously  underestimate 
the  heat  input  to  the  boundary  layer. 

At  the  present  time  the  Navy  ice  prediction  models  are 
only  intended  to  serve  as  an  aid  to  hand  analysis.  The  wind 
drift  model  used  by  the  Navy  ice  forecasting  group  at  the 
Fleet  Weather  Facility,  Suitland,  Maryland,  gives  ice  speed 
as  a linear  function  of  surface  wind  speed  [Skiles,  1968]. 
Moving  to  the  right  of  the  wind,  the  drift  angle  is  given  as 
an  inverse  exponential  function  of  wind  speed.  Under  de- 
velopment is  a model  to  provide  daily  wind  drift  after 
Shuleikin  [1953]  and  growth  after  Zubov  [1945]  at  62  coastal 
weather  reporting  stations  [Gerson  and  Simpson,  1976,  and 
Gerson,  1975].  Barnett  [1976]  has  presented  a scheme  for 
estimating  the  general  ice  conditions  during  the  shipping 
season  of  the  North  Slope  sea  route  using  the  height  of  the 
1000  mb  surface  at  two  geographic  locations  (52N,  100E,  and 
70N , 140E) . 
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All  these  models  are  conceptually  extremely  simple  in 
the  light  of  our  present  understanding  of  the  complexity  of 
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sea-ice  dynamics  and  they  will  probably  not  be  capable  of 
meeting  the  expanding  needs  of  the  Navy  for  arctic  ice  fore- 
casting. Several  more  complete  dynamical  models  of  the  arctic 
ice  pack  have  been  formulated  [Campbell,  1965;  Coon,  et  al, 
1977;  Hibler,  1977].  Therefore,  a new  model  is  in  order, 
but  first  an  assessment  of  present  Navy  analysis  and  forecast 
fields  in  the  Arctic  is  needed  to  determine  if  these  are 
capable  of  meeting  the  needs  of  the  more  complex  models. 


B.  DESCRIPTION  OF  THE  ARCTIC  ICE  PACK 

The  Arctic  ice  pack  forms  a lid  over  most  of  the  North 
Polar  Sea,  constantly  changing  its  extent  of  coverage, 
patterns  of  motion  and  thickness  distribution.  The  seasonal 
range  of  ice  coverage  is  illustrated  in  Figure  1.  Of  note 
is  the  extreme  variability  of  the  position  of  the  ice  edge, 
resulting  in  a 20%  change  in  the  areal  coverage,  receding  in 
summer  far  enough  to  allow  a shipping  season  along  the  north 
coast  of  Siberia  and  the  North  Slope  of  Alaska. 

The  summer  season  brings  melting  to  the  central  pack, 
reducing  the  average  thickness  by  about  one  meter  and  pro- 
ducing large  open  water  areas  [Orvig,  1970].  These  open 
water  areas  effectively  reduce  the  large  scale  strength  of 
the  ice  pack  to  zero  and  allow  the  floes  to  move  freely. 

While  the  day-to-day  movement  of  individual  ice  floes 
resembles  a random  walk,  the  long  term  general  circulation  of 
the  ice  pack  has  a regular  pattern.  Two  dominant  features 
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are  a trans-polar  drift  from  the  East  Siberian  Sea,  across 
the  North  Pole,  and  out  of  the  Arctic  into  the  East  Greenland 
Current,  and  an  anti-cyclonic  gyre  in  the  Beaufort  Sea.  These 
large  scale  climatic  motions  resemble  the  atmospheric  climatic 
circulation  patterns  in  the  same  way  that  the  lower  latitude 
ocean  circulations  resemble  climatic  wind  patterns  [Campbell, 

1965]  . 

In  contrast  to  the  relatively  free  movement  of  the  ice 
in  the  central  pack,  the  coastal  ice  remains  attached  to 
land  in  winter,  preventing  longshore  movement.  The  transition 
from  zero  velocity  of  the  land  fast  ice  to  a velocity  of  a 
few  centimeters/sec  in  the  central  pack  occurs  in  the  shear 
zone,  normally  found  a few  kilometers  offshore  where  sharpe 
ice  velocity  gradients  occur. 

The  central  pack  has  an  average  thickness  of  about  three 
meters  and  it  is  made  up  of  individual  floes  with  dimensions 
of  a few  meters  to  several  hundred  meters  [Orvig,  1970] . In 
winter  these  floes  are  separated  by  ridges  up  to  several 
meters  high  and  interlaced  with  leads  of  open  water  up  to 
0.5  km  wide  and  several  kilometers  long,  leads  which  quickly 
refreeze  to  form  large  smooth  areas  of  young  ice  which  even- 
tually form  ridges  and  new  floes.  These  ridges  and  leads  are 
formed  by  the  relative  motion  of  the  ice  whose  principal 
driving  force  is  the  wind.  The  nature  of  ridging  is  extremely 
complex.  Reduced  to  a simple  description,  sea-ice  ridges  are 
rubble  piles  formed  by  collisions  and  shear  interactions 
between  adjacent  floes.  The  ridges  above  level  ice  are 
accompanied  by  keels  below  the  surface,  keels  which  may 
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extend  many  times  as  deep  as  the  ridge  is  high.  Since  on  a 
large  scale  the  ice  is  densely  fractured  it  has  little  tensile 
strength  and  readily  opens  when  subjected  to  a stretching  or 
divergent  strain  [Pritchard  and  Colony,  1976]. 


II.  HISTORICAL  BACKGROUND 


k. 


The  modeling  of  ice  dynamics  can  be  divided  into  three 
categories;  empirical  modeling,  small-scale  analytical  model- 
ing and  synoptic-scale  analytic  modeling.  The  empirical 
modeling  began  with  the  observations  by  Nansen  during  the 
drift  of  the  Fram  across  the  Arctic  Ocean  (1893-96)  , who 

observed  that  ice  drift  speed  was  about  2%  of  the  wind  speed 

l j 

and  generally  to  the  right  of  the  wind  direction  at  an  average 

angle  of  28°.  Sverdrup,  drifting  with  the  ice  pack  in  the 
Maud,  (1919-25),  found  a wind  factor  for  ice  of  0.014  to 
0.024  and  a drift  angle  from  18°  to  40°,  the  latter  when  the 
ice  was  loose  in  summer. 

Zubov  [1945]  further  developed  the  empirical  relationship 
of  wind  and  ice  drift.  His  analysis  of  the  drift  of  the 
Sedov  yielded  a wind  factor  of  0.015  and  drift  angle  of  29°, 
which  compare  favorably  with  Nansen.  Zubov  further  postulated 
that  ice  drift  follows  the  sea-level  isobars,  i.e.,  it  is 
parallel  to  the  sea  level  geostrophic  wind  because  the  counter- 
clockwise turning  of  the  wind  through  the  atmospheric  friction 
layer  is  approximately  equivalent  to  the  ice  drift  angle  to 
the  right  of  the  surface  wind. 

Early  work  in  analytic  modeling  of  small-scale  ice 
dynamics  was  done  by  Shuleikin  [1953]  who  derived  the  force 

| 

balance  and  motion  of  the  ice  due  to  the  wind  stress.  However, 
his  theories  were  limited  to  10/10  ice  coverage.  Fukutomi 
» [1948]  treated  ice  drift  by  considering  the  effects  of  ice 
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concentration,  roughnes::  and  ocean  currents  in  an  analytical 
model . 


Modern  high  speed  computers  have  made  it  practical  to 
develop  synoptic-scale  ice  drift  models  which  span  a broad 
range  in  complexity.  Simplest  are  the  wind  drift  models 
using  empirical  relationships.  An  example  of  the  use  of  wind 
drift  relationships  is  the  numerical  model  developed  by  Knodle 
[1964],  which  combines  Shuleikin's  and  Fukutomi ' s theories, 
yielding  a wind  drift  field  which  is  a function  of  ice  con- 
centration for  Baffin  Bay.  The  most  complex  are  analytical 
models  which  consider  the  interactions  within  the  pack  ice. 
Because  of  the  internal  stresses  transmitted  through  the  ice 
pack  a solution  which  neglects  the  interaction  of  the  ice 
cannot  completely  describe  the  dynamics. 

The  first  models  which  attempted  to  include  the  inter- 
action of  the  ice  assumed  it  to  have  a viscous  nature.  In 
formulating  his  ice  circulation  model  Campbell  [1965]  considers 
individual  ice  floes  as  "...  fluid  elements  with  an  average 
area  which  is  small  compared  with  the  ocean  area,  [thus]  the 
ice  sheet  may  be  viewed  as  a film  of  highly  viscous  fluid 
suspended  between  two  less  viscous  fluids,  air  and  water." 

While  Campbell  and  several  other  authors  have  produced  viscous 
models  which  describe  the  general  circulation  of  ice  in  the 
Arctic,  all  have  limitations,  the  most  notable  is  their  in- 
ability to  model  the  coastal  shear  zone. 

Hibler  [1977]  has  developed  a viscous  model  which  simu- 
lates large  scale  (100  km)  sea-ice  dynamics;  the  model  has 
been  run  for  simulations  8 years  in  length.  Since  this  model 
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is  not  hampered  by  numerical  stability  problems  it  is  able 
to  utilize  large  time  steps  (1  day)  and  thus  perform  long- 
term simulations,  while  still  retaining  the  ability  to  be  run 
in  smaller  time  and  space  scales.  This  model  incorporates 
an  ice  thickness  evaluation  which  combines  mean  thickness 
and  compactness.  Despite  his  model's  success  Hibler  states, 
"...  particular  elastic-plastic  rheologies  may  well  be  more 
useful  in  reproducing  many  detailed  small  scale  effects." 

One  of  the  most  recent  models,  which  takes  a departure 
from  viewing  the  ice  as  viscous  in  nature,  is  that  formulated 
by  the  Arctic  Ice  Dynamics  Joint  Experiment  (AIDJEX)  program. 
It  models  the  ice  as  an  elastic-plastic  material  whose  yield 
strength  depends  on  the  amount  of  thin  ice  present.  A major 
difference  between  the  plastic  and  fluid  models  is  that  the 
plastic  material  responds  independently  of  the  rate  of  the 
deformation  [Coon,  Hall  and  Pritchard,  1977].  In  simulations 
covering  periods  of  about  10  days  the  AIDJEX  model  has  pro- 
duced good  results  in  describing  both  large  and  small  scale 
features  of  ice  dynamics  and  depicting  the  coastal  shear 
zone  [Coon,  et  al,  1977].  Because  of  the  potential  success 
of  the  AIDJEX  model  in  addressing  these  important  Navy  re- 
quirements for  ice  forecasting,  it  was  selected  as  a means 
of  evaluating  the  Navy's  atmospheric  analysis  and  prognostic 
fields  for  use  in  an  ice  dynamics  model. 
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III.  DESCRIPTION  OF  THE  AIDJEX  MODEL 


The  AIDJEX  ice  model  consists  of  three  major  parts:  the 
momentum  equation,  a set  of  constitutive  laws  which  describe 
the  stress-strain  behavior  of  the  ice  in  the  horizontal  plane, 
and  the  ice  thickness  distribution  [Coon,  et  al,  1976]. 

The  first  element  of  the  model  is  the  momentum  equation 
which  states  that  the  acceleration  of  an  element  of  ice  of 
mass  (m)  is  related  to  the  air  stress,  (r  ),  the  water  stress 

cl 

(^  ) , the  coriolis  force,  and  the  divergence  of  internal 
*+  = 

stresses  ( 7 • o ).  The  momentum  equation  is  written: 


mv  = r + T 
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= ice  velocity 
= coriolis  force 
= coriolis  parameter 
= sea  surface  tilt  force 
= geostrophic  current 


Air  stress  is  computed  as  a quadratic  function  of  the 
geostrophic  wind  [Pritchard,  et  al,  1976]: 
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(2) 


p=  = air  density 
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CD  = geostrophic  drag  coefficient 


where 


B 

- a 


U = geostrophic  wind  velocity 

cosot  -sin  ot 

sinot  cos  ot 

ct  = counterclockwise  drag  angle 

Water  stress  is  represented  by  a quadratic  drag  law 

similar  to  that  used  in  the  atmospheric  boundary  layer.  The 

stress  on  the  lower  surface  of  ice  (r  ) is  a function  of  the 

w 


ice  velocity  (v)  relative  to  v : 
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where 


T = PC  I v - v |B(v-v) 
w w w 1 g 1 - g 


P = water  density 
w 1 


(3) 


C = water  drag  coefficient  = 5.5  x 10 
w 3 
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B = 


cos(tt-£)  -sin  (7T-/S) 
sin  (7 r-/J)  cos  (ir-fi) 

P = drag  angle  = 23° 

Units  are  in  the  CGS  system.  The  geostrophic  flow  is  taken 
to  be  given  by  the  long  term  observed  values  [Pritchard, 
et  al,  1976] . 

The  stress  deformation  constitutive  law  has  been  chosen 
as  an  elastic-plastic  response.  A plastic  model  possesses 
many  of  the  properties  thought  to  be  desirable  in  describing 
the  response  of  ice  to  wind  stress.  Specifically,  it  is  able 
to  simulate  rapid  variations  which  approximate  discontinuous 
behavior  in  nearshore  regions  and  simulate  flow  resembling  a 
viscous  material  when  a large  enough  load  is  applied.  The 
constitutive  law  relates  internal  ice  stress  ( d ) to  the 


20 


•-V 


large  scale  strain  in  the  ice  sheet.  The  ice  is  a stiff 
elastic  material  allowing  elastic  deformation  of  only  about 

0. 1%  before  reaching  a yield  strength  in  compression  and 
shear.  Stresses  are  restricted  to  have  no  positive  principal 
values  because  the  material  is  densely  fractured  and  will  not 
support  stress  in  any  opening  mode  [Colony,  1976]. 

The  thickness  distribution  measures  the  percentage  of 
area  covered  by  ice  in  each  differential  range  of  thickness 
within  each  element  of  the  grid.  A typical  example  is  shown 
by  the  solid  line  in  Figure  2.  The  thickness  distribution 
determines  the  mass  per  unit  area  of  the  element  as  well  as 
its  mechanical  properties.  The  thickness  distribution  is 
affected  by  both  thermodynamic  and  mechanical  processes. 

For  example,  the  melting  of  a few  centimeters  off  the  top 
surface  would  shift  the  thickness  distribution  to  the  left, 

1. e.,  to  the  dashed  line  in  Figure  2.  During  the  winter,  all 
but  the  very  thickest  ice  grows  thicker,  shifting  the  curve 
to  the  right.  When  the  element  converges  slightly,  the  thin 
ice  is  ridged  into  thick  ice,  causing  the  shape  of  the  distri- 
bution to  change  as  some  thin  ice  disappears  and  a smaller 
increment  of  new  thick  ice  is  produced,  as  given  by  the 
dotted  line  in  Figure  2.  If  the  element  diverges,  open  water 
forming  between  the  floes  would  be  shown  in  the  distribution 
by  a delta  function  at  zero  thickness.  As  time  passes  this 
open  water  would  freeze  into  thin  ice,  moving  the  delta  function 
to  the  right  [Coon,  et  al,  1976].  In  the  AIDJEX  model  the 
growth  rate  is  defined  as  the  monthly  climatological  average 
[Pritchard  and  Colony,  1976]. 
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One  of  the  desirable  properties  of  the  AIDJEX  model  is 
the  potential  to  express  yield  strength  as  a function  of  the 
thickness  distribution.  However,  thus  far,  a satisfactory 
relationship  has  not  been  found  and  the  yield  strength  is 
varied  as  a unique  input  parameter.  The  redistribution 
function  and  the  energetics  argument  that  enables  strength 
to  be  determined  from  the  thickness  distribution  will  be 
reformulated  as  further  work  provides  direction  [Pritchard, 
et  al,  1977] . 

The  computational  grid  with  manned  camps  and  buoy  loca- 
tions for  15  May  1975  is  shown  in  Figure  3.  The  grid  was 
assembled  according  to  the  following  rules:  (1)  The  boundary 
may  be  a coastline  or  a dynamic  boundary  determined  by  the 
buoy  positions,  etc.  Positions  of  the  data  buoys  and  the 
manned  drifting  stations  were  found  using  the  Navy  Navigational 
Satellite  System  and  the  buoy  ring  motion  determines  the 
boundary  velocity  conditions.  This  is  the  only  boundary  con- 
dition used  in  the  model.  (2)  All  interior  cells  must  be 
quadrilaterals  and  all  cells  next  to  the  boundary  must  be 
quadrilaterals  or  triangles.  To  minimize  computer  time, 
linear  dimensions  of  all  cells  were  made  approximately  equal. 
The  stability  of  the  difference  scheme  depends  on  the  Courant 

, At  s'  ^ 

restriction,  C -j  t which  means  that  the  time  step  is  pro- 
portional to  the  cell  dimension,  specifically  the  dimension 
of  the  smallest  cell.  The  factor,  C,  is  a function  of  ice 
strength;  as  ice  strength  increases  the  time  step  must  de- 
crease . 
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Initial  conditions  for  a model  run  are  contained  on  a 


1 

computer  file,  magnetic  tape  or  disk,  which  must  be  generated 
from  a history  of  ice  motion  derived  from  buoy  and  manned 
camp  positions.  This  file,  called  a DUMP  file,  contains 
positions  of  the  grid  points,  initial  velocities  of  grid  cells, 
stress  state,  and  ice  thickness.  Periodically  during  the 
model  run  this  file  is  regenerated,  providing  a set  of  initial 
conditions  to  restart  the  run.  Thus  the  model  can  be  stopped 
and  restarted  several  times  during  a lengthy  run  to  allow 
checking  of  the  solution  or  to  rerun  the  problem  and  look  at 
some  event  on  a finer  time  scale.  For  example,  if  the  solution 
were  being  saved  for  plotting  every  four  hours,  it  might  be 
desirable  to  return  and  look  at  the  solution  every  twenty 
minutes  for  some  four-hour  period.  The  DUMP  file  also  serves 
as  the  data  source  for  plotting  two-dimensional  fields. 

One  additional  output  from  the  model  is  a time  history 
of  data  from  20  selected  nodes.  Saving  the  solution  at  only 
a fraction  of  the  nodes  enables  one  to  plot  the  solution 
densely  in  time,  enabling  high  frequencies  to  be  seen  without 
using  excessive  amounts  of  storage  [Don  Thomas,  AIDJEX, 

} 

personal  communication] . 

The  parts  of  the  model,  i.e.,  the  momentum  equation  and 
the  constitutive  laws  are  solved  simultaneously.  Detailed 
equations  and  methods  of  solution  are  given  by  Pritchard 
and  Colony  [1976] . 
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IV.  EVALUATION  OF  FNWC  SEA-LEVEL  PRESSURE  FIELD  AND 
GEOSTROPHIC  WINDS 

Fleet  Numerical  Weather  Central  (FNWC) , located  in 
Monterey,  California,  is  responsible  for  providing  all  numer- 
ical environmental  forecasts  in  support  of  Naval  operations. 
FNWC  is  linked  to  the  National  Meteorological  Center  (NMC) , 
the  Air  Force  Automated  Weather  Network  (AWN)  and  in  addition 
receives  observations  from  Navy  units  throughout  the  world 
[U.  S.  Naval  Weather  Service,  19751.  The  numerical  products 
produced  at  FNWC  include  sea-level  pressure,  upper-air  analyses 
and  prognoses,  marine  layer  winds,  sea-surface  temperature, 
and  many  others  produced  on  a polar  stereographic  projection 
from  data  analyzed  on  a rectangular  grid.  Most  fields  are 
resolved  on  a 63  x 63  point  Northern  Hemispheric  grid  but 
other  grids  are  used  for  some  products.  Atmospheric  prognoses 
(forecasts)  are  made  twice  a day  at  0000GMT  and  1200GMT  for 
each  12-hour  period  out  to  72  hours.  Many  computational 
algorithms  are  used  to  generate  the  numerical  products  but 
one  of  the  most  important  is  that  which  blends  available 
information  to  formulate  the  best  possible  sea-level  pressure 
field.  This  is  known  as  the  Fields  by  Information  Blending 
(FIB)  technique. 

A.  DESCRIPTION  OF  THE  FIELDS  BY  INFORMATION  BLENDING  (FIB) 

TECHNIQUE 

At  FNWC,  sea-level  pressure  is  analyzed  to  a 125  x 125 
polar  stereographic  projection  hemispheric  grid  (grid  spacing 
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190  km  at  60°N)  using  the  FIB  technique  [Holl  and  Mendenhall, 
1972].  While  analyzed  on  a 125  x 125  grid  the  resultant  field 
is  archived  on  a 63  x 63  grid  which  has  a 381-km  spacing  at 
60°N. 

The  FIB  generation  of  the  sea-level  pressure  analysis  is 
illustrated  in  Figure  4.  The  basic  steps  in  the  FIB  analysis 
technique  are : 

1.  The  6-hour  old  analysis  is  extrapolated  to  the 
present  analysis  time,  assembled  at  each  grid  point 
and  combined  with  the  6-  or  12-hour  old  prognosis 
to  produce  a first  guess  field. 

2.  New  wind  and  pressure  reports  are  assembled,  re- 
jecting errors  and  weighting  individual  values 
according  to  their  reliability. 

3.  Available  information  fields  are  blended  into  a 
modified  first  guess  field. 

4.  Fields  of  sea-level  pressure,  sea-level  pressure 
gradient  and  the  Laplacian  of  the  pressure  field 
are  analyzed  and  blended  to  obtain  the  best  fit 
field. 

5.  New  reliability  weights  are  determined. 

6.  A new  gross  error  check  and  reevaluation  of  data 
is  made. 

7.  Assembly,  blending  and  resultant  weight  operations 
are  recycled. 

8.  The  final  best-fit  pressure  field  is  output. 

Since  the  resulting  pressure  field  is  a blend  of  avail- 
able information,  in  the  absence  of  observations  the  analysis 
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will  be  largely  based  on  the  previous  prognosis.  However, 
the  quality  of  the  prognosis  is  dependent  upon  the  accuracy 
of  the  previous  analysis.  Thus,  in  data  sparse  areas  such 
as  the  Arctic,  the  analysis  is  prone  to  greater  error  than 
a data  rich  area. 

B.  SEA- LEVEL  PRESSURE  EVALUATION  METHOD 

Evaluation  of  FNWC's  atmospheric  sea-level  pressure 
analysis  was  performed  by  using  AIDJEX  sea-level  pressure 
observations  as  a standard.  Twenty-day  files  of  sea-level 
pressure  observations  and  geostrophic  wind  analyses  were 
obtained  from  AIDJEX  on  magnetic  tape.  The  period  from  1 May 
through  19  July  1975  (AIDJEX  days  121-200,  four  20-day  inter- 
vals) was  selected  for  detailed  evaluation.'*'  AIDJEX  data 
represent  high  density,  high  quality  observations  from  several 
manned  camps  and  remote  buoy  stations  whose  positions  on 
18  May  1975  are  shown  in  Figure  5.  The  FNWC  sea-level  pres- 
sure analysis  for  this  time  is  shown  in  Figure  6. 

FNWC's  6-hourly  analyses  and  12-hourly  prognoses  of  sea- 
level  pressure  fields  were  interpolated  to  each  AIDJEX  ob- 
servation station.  The  RMS  pressure  difference  (where  pres- 
sure difference  is  ^ p=pfnwC-PAIDJEX*  at  AIDJEX  stations 
in  each  20-day  interval  was  calculated  for  the  analysis,  24-, 
48-,  and  72-hour  prognoses,  and  the  results  are  summarized 
in  Figure  7.  The  RMS  pressure  difference  (a  measure  of  the 


^"AIDJEX  days  are  Julian  days  measured  from  1 Jan  1975 
and  continue  to  be  consecutively  numbered  into  1976.  Thus 
2 Jan  1976  is  AIDJEX  day  366. 


quality  of  the  product)  increases  from  2.3  mb  in  the  first 
interval  (A,  1-20  May)  to  3.4  mb  in  the  last  interval  (D, 

30  June  to  19  July) . For  the  prognostic  times  the  error 
increases  by  about  0.6  mb  for  the  first  24-hour  period  and 
thereafter  by  about  2 mb  for  each  24-hour  projection.  For 
some  reason  the  second  and  third  time  intervals  (B,  21  May  to 
9 June  and  C,  10-29  June)  show  an  increase  in  quality  over 
the  first  interval.  However,  the  quality  is  always  poorest 
in  the  final  interval  which  is  the  summer  period. 

An  explanation  for  the  larger  RMS  4 P values  in  the 
fourth  time  interval  is  found  in  Figure  8 which  shows  a time 
series  of  atmospheric  sea-level  pressure  at  the  AIDJEX  main 
camp  over  the  entire  90-day  period.  The  greatest  variability 
of  the  pressure  is  in  the  fourth  period  implying  an  increased 
synoptic  activity  level.  This  increased  activity  combined 
with  lack  of  observations  probably  accounts  for  the  decreased 
quality. 

An  attempt  was  made  to  determine  if  and  when  observations 
from  AIDJEX  were  included  in  the  FNWC  analysis.  A search  was 
made  of  FNWC  archives  for  any  reports  submitted  by  AIDJEX. 

A total  of  95  AIDJEX  reports  were  found  in  the  archives,  re- 
ported only  from  the  main  AIDJEX  camp.  However,  due  to  the 
nature  of  the  archiving  system  late  reports  (up  to  24  hours) 
are  archived  and  it  is  not  possible  to  determine  if  any  par- 
ticular report  was  included  in  the  analysis.  RMS  4P  for 
analyses  corresponding  to  archived  reports  is  shown  by  the 
line  marked  R on  Figure  7.  A careful  examination  of  Figure 
3 reveals  the  effect  of  the  inclusion  of  an  AIDJEX  pressure 


observation  in  the  FNWC  analysis.  Time  of  archived  AIDJEX 
observations  are  indicated  by  "tick"  marks  near  the  lower 
edge  of  the  graph.  The  circled  areas  illustrate  a few  times 
corresponding  to  archived  observations  when  the  analysis  has 
made  a sharp  change  in  trend  to  become  coincident  with  the 
value  of  the  observations.  It  can  be  inferred  that  the  AIDJEX 
reports  did  reduce  errors  in  the  FNWC  analysis  by  0.7  mb  on 
the  average  and  that  any  observations  from  the  central  Arctic 
will  probably  result  in  an  improved  analysis. 

Walsh  [1977]  has  calculated  the  RMS  height  errors  for 
the  1000  mb  height  forecasts  of  the  NMC  primitive  equation 
model  at  16  grid  points  covering  roughly  the  same  area  as 
the  AIDJEX  observation  stations.  His  technique  involved 
finding  the  height  error  of  the  24-,  48-,  and  72-hour  fore- 
casts using  the  verifying  analysis  as  a standard.  When  con- 
verted to  millibars  (8  m = 1 mb)  these  errors  range  from 
2-3  mb  at  24  hours  to  5-6  mb  at  72  hours.  This  technique  of 
using  the  verifying  analysis  vice  observations  as  a standard 
was  applied  to  the  FNWC  data,  again  interpolated  to  the  same 
AIDJEX  station  locations.  The  RMS  A P increases  from  2.5  mb 
at  24  hours  to  6.5  mb  at  72  hours.  These  are  comparable  to 
the  results  found  by  Walsh.  However,  they  show  a marked 
improvement  over  the  use  of  observations  as  a standard,  the 
improvement  ranging  from  1.0  at  24  hours  to  .3  mb  at  72  hours. 
This  suggests  that  the  use  of  the  verifying  analysis  as  a 
standard  for  forecast  evaluation  yields  optimistic  results 
primarily  because  the  analysis  is  largely  based  on  the  pre- 
vious prognosis  in  the  Arctic  area. 
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Sufficient  pressure  data  were  available  from  AIDJEX  to 
evaluate  FNWC  analyses  from  1 May  1975  through  24  April  1976. 
The  results  of  the  evaluation  are  shown  in  Figure  9.  d P 
began  at  about  0.3  mb  (higher  pressure  than  actually  existed) 
and  continued  increasing  into  the  summer  until  the  end  of 
August  when  the  value  reached  2.7  mb.  It  decreased  as  fall 
approached  becoming  negative  in  mid-October,  and  generally 
remained  negative  throughout  the  winter,  becoming  positive 
the  following  March. 

The  RMS  dP  remained  high  throughout  the  summer,  de- 
creased slightly  in  the  fall  and  remained  lower  than  the  sum- 
mer values  throughout  the  winter.  The  sharp  increase  in 
March  can  be  attributed  to  increased  synoptic  activity.  An 
evaluation  of  whether  these  fluctuations  of  RMS  dP  and  dP 
follow  an  annual  cycle  will  become  possible  when  a longer 
record  of  complete  observations  in  the  Arctic  become  avail- 
able . 


C.  METHOD  OF  EVALUATION  OF  GEOSTROPHIC  WINDS 

An  evaluation  of  FNWC  geostrophic  winds  was  performed  to 
determine  their  suitability  for  air  stress  computation.  In 
the  AIDJEX  model  and  probably  any  other  operational  model  the 
wind  stress  is  determined  from  the  geostrophic  wind.  There- 
fore, the  quality  of  the  sea-level  pressure  field  is  an  im- 
portant consideration.  The  AIDJEX  sea-level  pressure  field 
is  expressed  by  a least  squares  polynomial  fit  to  all  data 
points  [AIDJEX,  1976].  However,  geostrophic  winds  are  de- 
pendent on  the  pressure  gradient  rather  than  the  actual 


pressure  values  at  any  point.  Geos trophic  winds,  as  obtained 
from  AIDJEX,  are  found  by  differentiating  the  polynomial  ex- 
pression along  grid  directions  at  each  observation  station. 

To  obtain  FNWC  geostrophic  winds,  a fourth  order  finite 
differencing  technique  was  applied  to  the  FNWC  grid  point 
sea- level  pressure  value  to  determine  the  pressure  gradient 
components.  FNWC  geostrophic  wind  components  were  then  inter- 
polated to  AIDJEX  observation  stations  and  compared  with 
AIDJEX  geostrophic  winds.  The  RMS  wind  speed  and  direction 
differences  were  calculated  and  the  results  are  summarized 
in  Figure  10  where  wind  speed  differences  have  been  normalized 
by  the  AIDJEX  value.  Aggregate  results  of  all  four  evaluation 
periods  are  marked  by  *G  which  shows  increased  differences 
as  the  forecast  period  increases.  In  particular  the  direction 
differences  increase  from  about  50°  to  75°  from  analysis  to 
the  72-hour  forecast,  respectively.  The  fractional  speed 
errors  show  an  increase  from  0.65  to  1.35  for  the  same  periods. 

If  only  wind  speeds  in  excess  of  5 m/sec  are  considered, 
the  RMS  differences  decrease  significantly,  and  only  a slight 
increase  from  0.4  to  0.5  is  found  as  the  forecast  period 
lengthens.  These  values  are  of  particular  interest  since  the 
air  stress  is  taken  as  proportional  to  the  square  of  wind 
speed  and  the  higher  wind  speeds  will  be  of  greater  importance 
in  determining  ice  drift.  There  is  a relative  minimum  in  both 
direction  and  speed  differences  in  the  21  May  to  9 June  time 
period  which  corresponds  to  the  reduced  RMS  pressure  error 
in  this  time  frame.  However,  for  some  reason  the  direction 
and  speed  difference  for  period  C do  not  similarly  correspond 
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with  the  relatively  low  RMS  pressure  errors  for  that  time 
frame. 

D.  SIMULATION  RESULTS 

To  more  completely  evaluate  the  usefulness  of  the  FNWC 
sea-level  pressure  field  in  driving  an  ice  dynamics  model, 
a simulation  using  the  AIDJEX  model  was  made.  The  AIDJEX 
model  was  programmed  on  the  FNWC  6500  computer  to  simulate 
ice  movement  during  the  period  0300GMT  15  May  to  1200GMT 
25  May  1975.  Air  stress  derived  from  the  FNWC  sea-level 
pressure  analysis  field  was  used  in  the  simulation  while  all 
other  input  parameters  and  boundary  conditions  were  specified 
as  those  used  by  AIDJEX.  The  air  stress  field  was  computed 
using  equation  (2)  and  assembled  in  a form  compatible  with 
the  AIDJEX  model  using  the  following  values  of  air  density 
( Pa) , drag  coefficient  {CQ)  and  turning  angle  ( a ) : 

P = 0.00135  gm/cm3 

cL 

CD  = 11  x 10~4 
a = 20° 

CD  and  a are  the  same  values  used  in  AIDJEX  model  simulations 
[Eric  Levitt,  AIDJEX,  personal  communication] . 

Figures  11  and  12  show  time  series  plots  of  both  FNWC 
and  AIDJEX  computer  air  stress  magnitude  and  direction  at  two 
grid  points.  It  is  clear  that  the  FNWC  air  stress  is  smaller 
than  the  AIDJEX  air  stress,  in  particular  for  the  important 
time  periods  days  135-136  and  137-140,  when  there  is  a high 
air  stress.  At  node  7,7,  which  corresponds  to  the  main 
AIDJEX  camp,  during  day  135  the  FNWC  air  stress  is  20%  of 


AIDJEX  values.  From  day  137  to  140  the  FNWC  air  stress  is 
in  general  much  less  than  AIDJEX.  The  values  at  node  3,7 
follow  a similar  pattern,  the  FNWC  stress  being  a small 
fraction  of  the  AIDJEX  value  for  the  first  day,  and  the  FNWC 
stress  generally  retaining  the  trends  of  the  AIDJEX  stress, 
but  at  reduced  levels.  In  addition  to  the  magnitude  differ- 
ences there  are  direction  differences  of  up  to  150°,  however, 
these  are  generally  limited  to  the  light  wind  conditions. 

The  direction  differences  for  the  critically  high  stress 
period  137-140  are  small,  i.e.,  generally  less  than  30°. 

Figures  13  and  14  show  vector  field  plots  of  air  stress 
at  1200 GMT  18  May  for  the  AIDJEX  and  FNWC  model  runs,  respec- 
tively. The  AIDJEX  plot  shows  considerable  spatial  variation, 
with  an  area  of  high  stress  over  the  southern  portion  of  the 
region.  The  FNWC  stress  field  plot  shows  less  spatial  vari- 
ation and  significantly  smaller  air  stress  over  the  north- 
eastern portion  of  the  region. 

The  trajectories  computed  using  FNWC  and  AIDJEX  stress 
fields  for  the  main  AIDJEX  camp  are  shown  in  Figure  15.  The 
total  displacements  of  the  FNWC  projected  movement  is  compar- 
able to  the  AIDJEX  movement  with  only  a 3 km  difference  in 
the  final  positions.  In  addition  to  the  comparable  total 
displacement,  the  small-scale  features  of  the  AIDJEX  movement 
are  retained  in  the  FNWC  movement.  Figure  15  also  includes 
the  actual  movement  of  the  AIDJEX  main  camp  as  calculated 
from  AIDJEX  measured  positions  [Thorndike  and  Cheung,  1977]. 

It  is  clear  that  neither  AIDJEX  nor  FNWC  followed  the  movement 
exactly.  During  the  first  several  days  the  FNWC  movement 
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more  closely  resembles  the  true  motion,  while  in  the  final 
positions  the  AIDJEX  simulation  is  closer,  namely  3 km  versus 
4.5  km. 

A more  detailed  analysis  is  possible  by  looking  at  the 
x-  and  y-position  time  series  for  the  main  camp.  Figure  16 
shows  the  x-  and  y-position  time  series  relative  to  the  AIDJEX 
coordinate  system  [Thorndike  and  Chueng,  1977].  It  is  evident 
that  the  error  in  the  y-position  was  accumulated  largely  in 

the  first  two  days.  Referring  back  to  Figure  11,  the  FNWC 

. 2 
air  stress  is  only  about  0.3  dynes/cm  compared  to  1.5  dynes/ 

2 

cm  for  AIDJEX.  This  low  FNWC  value  is  the  cause  of  most  of 
the  difference  in  the  final  position. 

The  relative  motion  of  two  or  more  stations  is  important 
for  those  processes  such  as  deformation  which  depend  on 
differential  motion.  The  trajectory  of  node  3,7,  located  in 
the  left  hand  portion  of  the  grid,  is  shown  in  Figure  17.  It 
has  a larger  scale  of  motion  than  the  main  camp  and  the  final 
error  is  in  the  opposite  direction,  being  about  3 km  in  the 
negative  y direction.  This  quantitatively  implies  a larger 
deformation  difference  between  the  two  stations  than  the 
absolute  motion  difference  of  either  single  station.  Node 
3,7  had  no  position  measurements  to  allow  checking  of  the 
results  of  the  simulations. 

Figure  18  shows  the  x-  and  y-position  time  series  for 
node  3,7.  Unlike  node  7,7,  whose  error  was  accumulated 
largely  in  the  first  day,  this  node  accumulated  its  error  in 
the  138-139  day  time  frame.  For  node  3,7  the  error  steadily 
increases  from  day  137-139.5,  this  time  period  corresponding 
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to  the  maximum  difference  in  air  stress  between  the  two  runs 
(Figure  11) . 

Figure  19  is  a vector  plot  of  the  ice  velocity  field 
for  the  AIDJEX  simulation  at  1200GMT  18  May.  Recalling  that 
the  boundary  velocity  is  specified,  it  is  evident  that  the 
discontinuities  in  the  velocity  field  along  the  southwest 
and  northeast  boundaries  represent  errors  in  the  calculated 
velocities.  Figure  20  is  a vector  plot  of  the  stress  diver- 
gence field  at  the  same  time  showing  relatively  large  stress 
divergence  forces  along  the  south  and  northeast  boundaries. 
The  turning  of  the  vector  along  the  south  boundary  implies 
that  the  stress  is  an  artifact  of  the  discontinuous  velocity 
field  and  probably  has  the  effect  of  directing  the  field 
velocity  into  agreement  with  the  boundary  velocity. 

Figure  21  is  a vector  plot  of  the  velocity  field  for 
the  FNWC  simulation  at  1200GMT  18  May.  The  velocity  dis- 
continuity along  the  northeast  boundary  is  similar  to  AIDJEX 
but  there  is  no  apparent  discontinuity  along  the  south  bound- 
ary. The  regular  pattern  of  the  stress  divergence  shown  in 
Figure  22  appears  to  be  a more  natural  occurrence  than  does 
the  AIDJEX  stress  divergence  field,  implying  that  at  least 
for  this  one  time  the  FNWC  simulation  is  a better  representa- 
tion than  the  AIDJEX  simulation  of  the  actual  ice  dynamics. 

The  relationship  of  the  two  computed  velocities  is  shown 
in  Figure  23.  The  only  large  differences  in  the  velocities 
are  for  days  135  and  137.  These  correspond  to  the  large 
differences  in  air  stress  values.  The  remaining  time  period 
shows  a high  correlation  between  the  two  data  sets,  an  even 
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higher  correlation  than  would  seem  to  be  supported  by  the 
significant  differences  in  the  air  stress  alone.  The  reasons 
for  this  high  correlation  can  be  attributed  to  the  influences 
of  the  boundary  velocity  condition.  The  ice  is  largely 
affected  by  internal  ice  stress  during  the  light  wind  con- 
ditions and  is  primarily  driven  by  the  boundary  velocity 
inputs.  This  is  amply  illustrated  in  Figures  24  through  27 
where  the  FNWC  and  AIDJEX  air  stress  fields  are  radically 
different  yet  the  velocity  fields  are  nearly  identical.  The 
model  is  being  run  by  the  boundary  at  this  time. 
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V.  CONCLUSIONS 


A comparison  of  FNWC  sea-level  pressure  analyses  and 
prognoses  with  observations  made  by  AIDJEX  during  the  spring 
and  summer  of  1975  shows: 

1.  A RMS  error  of  approximately  3 mb  for  the  analysis 
time  increasing  to  approximately  7 mb  at  72  hours. 

2.  The  RMS  error  is  highest  during  the  summer  months, 
increasing  from  2.3  mb  in  May  to  3.8  mb  in  July. 

Comparing  the  FNWC  sea-level  pressure  analyses  and  AIDJEX 
observations  for  the  360-day  period  from  1 May  1975  to  24 
April  1976  shows: 

1.  The  RMS  error  remains  relatively  high  throughout 
the  summer,  decreases  into  fall  and  winter  with  a 
sharp  rise  in  February. 

2.  The  average  4?  increases  to  a maximum  of  1.7  mb  in 
August,  becoming  negative  in  October  reaching  a 
minimum  of  -1.0  mb  in  November,  and  then  becomes 
positive  again  in  late  winter. 

Of  particular  note  is  the  fact  that  the  inclusion  of  just 
one  observation  for  the  AIDJEX  ocean  improves  the  pressure 
analysis  by  nearly  1 mb. 

A comparison  of  FNWC  generated  geostrophic  winds  with 
AIDJEX  generated  geostrophic  wind  reveals: 

1.  The  normalized  RMS  speed  difference  for  all  wind 
speeds  increases  from  0.65  to  1.35  from  analysis 
to  the  72-hour  prognosis. 
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2.  Considering  only  wind  speeds  greater  than  5 m/sec, 
the  normalized  RMS  speed  difference  is  lower  and 
increases  from  0.4  to  0.5  from  analysis  to  72-hour 
prognosis . 

The  movement  of  two  AIDJEX  grid  nodes  using  FNWC  sea- 
level  pressure  fields  to  derive  wind  stress  in  a ten-day  simu- 
lation closely  approximates  the  movement  calculated  by  AIDJEX. 
Final  position  errors  compared  to  measured  positions  of  the 
main  camp  were  4.5  km  for  the  FNWC  simulation  versus  3 km  for 
the  AIDJEX  simulation.  The  other  node  had  a difference  of 
about  3 km  between  AIDJEX  and  FNWC  simulations.  The  majority 
if  not  all  the  final  position  error  was  accumulated  during 
periods  of  high  air  stress  when  FNWC  air  stress  was  20%  to 
80%  of  the  AIDJEX  value. 

During  periods  of  low  wind  stress  the  AIDJEX  model  is 
largely  driven  by  the  boundary  velocity  conditions  as  evidenced 
by  the  high  correlation  of  AIDJEX  and  FNWC  calculated  veloci- 
ties during  these  periods. 

The  performance  of  the  AIDJEX  model  is  highly  dependent 
on  the  accurate  specification  of  the  boundary  velocities. 

Any  attempt  to  use  the  model  as  a predictive  tool  will  be 
limited  by  the  use  of  the  boundary  conditions.  The  present 
method  is  not  realistic  for  forecasting  and  must  be  replaced 
by  some  form  of  wind  driven  boundary.  A significant  degrada- 
tion of  the  performance  of  the  model  can  be  expected  to  result 
from  the  use  of  a wind  driven  boundary. 
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Another  limitation  of  the  model's  use  for  forecasting 
is  the  amounts  of  computer  time  required  for  a simulation. 

A 72-hour  simulation  required  about  1.5  hours  computer  time 
on  a CDC  6500  computer.  This  is  fully  one-fourth  the  time 
required  to  generate  a 72-hour  7-level  primitive  equation 
meteorological  prognosis  for  the  entire  Northern  Hemisphere. 

The  important  ice  thickness  distribution  was  not  investi- 
gated in  the  simulation  because  it  was  deleted  from  the  model 
for  reasons  relating  to  its  contribution  to  the  ice  strength 
calculations.  However,  it  may  be  included  later  if  one 
desires  it  for  other  than  ice  strength  considerations. 
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VI.  RECOMMENDATIONS 

It  is  recommended  that  the  Navy  continue  work  on  opera- 
tional ice  dynamics  models.  Specific  areas  of  investigation 
should  include  a formulation  of  boundary  conditions  which 
are  compatible  with  a forecasting  vice  hindcasting  scheme, 
preparation  of  a model  which  works  well  in  the  summer  ship- 
ping season,  and  incorporation  and  evaluation  of  an  ice 
thickness  distribution  technique  which  is  independent  of  ice 
strength  calculations. 

Furthermore,  every  effort  should  be  made  to  provide 
synoptic  sea-level  pressure  and  air  temperature  observations 
from  the  central  Arctic.  The  satellite  telemetry  buoys 
developed  and  used  in  conjunction  with  the  AIDJEX  field 
experiment  offer  one  of  the  simplest  and  most  direct  means 
of  improving  the  environmental  forecasts  in  the  Arctic.  With- 
out the  real  time  availability  of  these  observations,  good 
predictions  of  ice  conditions  encompassing  all  Navy  require- 
ments will  remain  an  impossibility. 
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APPENDIX  A FIGURES 


Figure  1.  Extent  of  sea  ice  in  the  Northern  Hemisphere.  [Sater, 
Ronhovde,  and  Van  Allen,  1971] 
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Figure  2.  A representative  thickness  distribution  undergoing  changes 
described  in  the  text.  [Coon,  et  al , 1976] 
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Figure  3.  The  AIDJEX  computational  grid  at  0300GMT  15  May  1975. 
Observation  stations  are  marked  as:  remote  buoys  (#  ),  manned 
camps  ( ■ ),  main  AIDJEX  camp  (node  7,7)  ( ■ BB) . Node  3,7  is 
marked  as  ( A ) • Irregular  configuration  of  the  grid  is  a result 
of  applying  the  rules  for  generating  the  grid  as  described  in  the 
text . 
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Figure  4.  Simplified  flow  diagram  of  FNWC's  Fields  by  Information 
Blending  (FIB)  Sea-Level  Pressure  analysis  technique. 


Figure  S.  AIDJEX  sea-level  pressure  analysis  for  1200GMT  18  May  1975. 
Values  are  sea- level  pressure  minus  1000  mb.  AIDJEX  observation 
stations  are  marked  by  £ ■ ) . Station  marked  B3  is  the  location  of 
AIDJEX  main  camp  (Big  Bear) . 


Figure  6.  FNWC  sea- level  pressure  analysis  for  1200GMT  18  May  1975. 
63x63  grid  points  are  marked  ( # ) . Station  marked  ( ■ BB)  is  location 
of  AIDJEX  main  camp.  Values  are  sea- level  pressure  minus  1000  mb. 
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Figure  7.  RMS  pressure  error  for  FNWC  sea-level  pressure  analysis 
and  prognoses  at  all  AIDJEX  stations.  Analysis  time  is  T=0,  T=24  is 
24  hour  prognosis,  T=48  is  48  hour  prognosis,  T=72  is  72  hour  prog- 
nosis. Dates  are  shown  by  A,  1 May  to  20  May;  B,  21  May  to  9 June; 

C,  10  June  to  29  June;  and  D,  30  June  to  19  July.  Points  marked  *G 
are  aggregate  errors  for  all  four  time  periods.  Line  marked  R is  RMS 
error  for  only  those  analysis  times  corresponding  to  FNWC-archived 
AIDJEX  observations. 
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Figure  9.  RMS  pressure  errors  and  mean  pressure  errors  for  FNWC 
analysis  at  all  AIDJEX  stations  for  the  period  1 May  1975  to  24  April 
1976.  Dates  indicated  are  AIDJEX  days.  Each  point  represents  20-day 
average  value.  The  oeginning  of  each  month  is  indicated.  RMS  error 
is  shown  by  a solid  line,  mean  error  by  a dashed  line. 
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Figure  10.  RJi  wind  speed  and  direction  difference  between  AIDJEX 
anal/ted  and  FNWC  analyzed  and  prognosticated  winds.  Lines  marked  5 are 
differences  for  only  those  wind  speeds  greater  than  5 m/sec.  Points 
marked  *G  are  aggregate  differences  for  all  four  time  periods.  Time 
periods  are  A,  1 May  to  20  May;  B,  21  May  to  9 June;  C,  10  June  to 
29  June;  and  D,  30  June  to  19  July. 


Figure  11a.  Computed  air  stress  magnitude  (TAUA-MAG)  from  AIDJEX 
field  (solid  line)  and  FNWC  field  (dashed  line)  at  AIDJEX  main  camp 
for  days  15  to  25  May  1975.  Units  are  dynes/cm^. 
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Figure  12a.  Computed  air  stress  magnitude  (TAUA-MAG)  from  AIDJEX 
field  (solid  line)  and  FNWC  field  (dashed  line)  at  node  3,7  for  days 
15  to  25  May  1975.  Units  are  dynes/cm2. 


Figure  12b.  Computed  air  stress  direction  (TAUA-DIR)  from  AIDJEX 
field  (solid  line)  and  FNWC  field  (dashed  line)  at  node  3,7  for  days 
15  to  25  May  1975.  Units  are  degrees  true. 
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Figure  13.  AIDJEX  computed  air  stress  field  at  1200GMT  18  May  1975. 

Scale  vector  represents  6 dynes/cnA. 
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Figure  14.  FNWC  computed  air  stress  field  at  1200GNTT  18  May  1975. 
Scale  vector  represents  6 dynes/cm^. 
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Figure  15.  Trajectories  of  AIDJEX  main  camp  from  15  to  25  May  1975 
for  AIDJEX  simulation  (solid  line)  and  FNWC  simulation  (dashed  line). 
( ■ ) Starting  position,  ( # ) ending  position,  ( • ) actual  measured 
position  every  six  hours  with  AIDJEX  dates  indicated.  Grid  spacing 
is  5 km/div. 
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Figure  16b.  Y - Position  time  series  of  AIDJEX  main  camp  as  computed 
by  AIDJEX  simulation  (solid  line)  and  FNWC  simulation  (dashed  line) . 
Ordinate  spacing  5 km/div. 


a.yeit 


X-POSITION  VS  Y-POSITION* CYCLE  0 TO  4981.  NODE 


Figure  l".  Trajectories  of  node  3,7  from  15  to  25  May  1975  for  AIDJEX 
simulation  (solid  line)  and  FNWC  simulation  (dashed  line).  Grid 
spacing  5 km/div. 
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VEL( N+l/2 ) CYCLE=1621  TIME=138/12  0/  0 


SCALE  30 . 00.-" 


75RUN1F/P*=4E7/HBAR=300/G=0  /INITIAL  STRESS 


P=.  67 


Figure  19.  AIDJEX  computed  velocity  field  at  1200GMT  18  May  19/5 

Scale  vector  represents  30  cm/sec. 


621  TIME=138/12  2/  0 


SCALE 


75RUN1F/P*=4E7/HBFIR=300/G=0  /INITIAL  STRESS  P=.67 


AIDJEX  computed  stress  divergence  field  at  1200GMT 
Scale  vector  represents  2 dynes/cm^. 


VELC N+l/2 ) CYCLE=1621  TIME=138/12  0/  0 


FNNC  RUN  135-2. USING  REVISED  FNWC  RIR  STRESS 


Figure  21.  FNWC  computed  velocity  field  at  1200GMT  18  May  1975. 
Scale  vector  represents  30  cm/sec. 
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VEI.-nZI  VS  TIME.  CYCLE  0 TO  4981,  NODE 


IZB-13A 


Figure  23b.  Velocity  direction  of  AIDJEX  main  camp  as  computed  by 
AIDJEX  simulation  (solid  line)  and  FNWC  simulation  (dashed  line). 
Direction  units  are  degrees  true. 


AIR  STRESS  CYCLE=3541 


'IME=142/12  0/  0 


75RUN 1F/P*=4E7/HBAR=300/G=0/ IN  I T I AL  STRESS  P- . 67P* 


Figure  24.  AIDJEX  computed  air  stress  field  at  1200GMT  22  May  1975. 
Scale  vector  represents  2 dynes/cm^. 


VEL(N+l/2)  CYCI_E=3541  TIME-142/12  0/  0 


75RUN1F/P*=4E7/HBRR=300/G=0/ INITIAL  STRESS  P-. 67P* 


Figure  25.  AIDJEX  computed  velocity  field  at  1200GMT  22  May  1975. 
Scale  vector  represents  10  cm/sec. 
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RIR  STRESS  CYCLE=3541  TIME=142/12  0/  0 


Figure  26.  FN’WC  computed  air  stress  field  at  1200GNTT  22  May  1975. 
Scale  vector  represents  2 dynes/cm^. 
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VEL( N+l/2)  CYCLE=3541  TIME 


SCALE  10.00. 


USING  REVISED  FNWC  AIR  STRESS 


FNWC  RUN135 


Figure  27.  FNWC  computed  velocity  field  at  1200GMT  22  May  1975 
Scale  vector  represents  10  cm/sec. 
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